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Abstract 

The pentaquark state recently discovered has been studied with three quark models which either 
fit the nucleon spin structure or the A^A^ scattering. A minimum 0"^ mass of 1620 MeV is obtained 
both for the state. The mixing of various color structure configurations, which would reduce 
the mass of the penta-quark state, should be taking into account in the calculation of penta-quark 
mass. 
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Eleven groups Q claimed recently that they found a penta-quark state, now called 9"*", 
with mass~1540 MeV, and width F < 25 MeV. This state is either identified from the decay 
product nK^ or pKs-, but up to now no experiment has identified both. In addition, the 
NA49 collaboration claimed that they found the anti-decuplet partner S of 0+ [2]. The 
HERA-Hl collaboration claimed that they found the charm penta-quark Be [3]. 

These measurements might be contaminated by the normal meson production due to the 
kinematical refiection 0,0]. The NA49 claim has been challenged by another CERN group 
based on S spectroscopy data with higher statistics |^. HERA-B p-nucleus reaction data 
has not found the G"*" [7J. BES J/\E' decay data analysis has not found the 0"*" either 
There are other groups that have not found the G"* 



The results of the reanalysis of the K^d and K^j^p scattering data are diverse lll |. 
W.R. Gibbs reanalyzed the K^d scattering data, taking into account the double scattering, 
and found a structure corresponding to a resonance with width of 0.9±0.2 MeV, which is 
either a 1.547±0.002 GeV \~ or a 1.559±0.003 GeV |^ state and \~ is favored In 
addition, a very tiny bump had appeared in 1973 CERN K^p pK^TT~^ inelastic scattering 
data 

Theoretical studies based on the chiral soliton quark model played an important role in 
triggering the G"*" searches In the chiral soliton quark model, the G^ is a member of the 
anti-decuplet rotational excitation following the well established octet and decuplet baryons 
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le QCD background of this model has been debated between different groups 
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Various quark models have been proposed to understand the G"*", mainly aimed at under- 
standing the parity, low mass and narrow width. In the naive quark model Q| the ground 
state should be an S-wave one and this means G"*" should be a negative parity state. The 
S-wave uudds configuration will have S'-wave KN components. However both the / = 0, 1 
KN S'-wave phase shifts are negative in the G+ energy region [2^ and this means that G"*" 
can not be an S'-wave KN resonance, except the new analysis of Gibbs On the other 
hand, since the / = KN P-wave Pqi phase shifts are positive, there might be resonance 
in this channel and this is consistent with the = predicted by the chiral soliton quark 



model. Bag model gave similar result 
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Various quark correlations have been proposed such as color Cooper pair and others 
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28l | aimed to reconcile the difference between the quark model predictions and the 
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properties of Q"*". It seems to be able to get an even p arity penta-quark ground state with 
small width but hard to get as low as 1540 MeV mass 29|. For example, Yu. A. Simonov 
did a quantitative calculation based on the Jaffe-Wilczek configuration |l7| by means of the 
effective Hamiltonian approach. The calculated G"*" mass is about 400 MeV higher than 



the observed 1540 MeV 



. Different quark interaction mechanisms such as quark-gluon, 



Goldstone boson exchange and instanton interaction have all been tried to understand the 



0+ 
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32| |. It is fair to say up to now there is no ab initio dynamical calculation to 
obtain a B"^ mass as low as 1540 MeV with the constraint that these model parameters fit 
the normal hadron spectroscopy 



QCD sum rule and lattice QCD both have been used to calculate the penta-quark, the 
ults are diverse. One Lattice QCD grou p r eported they have not observed either / = or 
= 1, = bound penta-quark state [35]. Two groups observed odd and one observed 



even parity state j^]. The pitfalls of these lattice QCD calculations have been discussed in 
j^^l • The is a resonance state and so its mass should be a complex number. It might be 
difficult to detect such a complex number in an Euclidean Monte Carlo calculation. QCD 
sum rule calculations are diverse too. Zhu and Sugiyama et al. favor negative parity ones 
while others favor positive parity 3si- Kondo et al. suggested to remove the reducible part 
from the correlation function and after do that they got positive parity penta-quark. 

Our group has done three quark model calculations. The first one is an application of 
the Fock space expansion model which we developed to explain the nucleon spin structure 
j^^. The naive quark model assumes that the baryon has a pure valence configuration. 
This is certainly an approximation. One expects there should be higher Fock components, 

B = aq^ + bq^qq H . (1) 

The nucleon spin structure discovered in polarized lepton-nucleon deep inelastic scattering 
shows that there are intrinsic non-perturbative sea quark components in nucleon and indeed 
the nucleon spin structure can be explained by a dynamical model of nucleon where the 
ground state has about 15% q'^qq component. This means that even the nucleon ground 
state might be a mixing of tri- and penta-quark components. In the mass calculation 
we assume it is a pure uudds five quark state but with channel coupling. Our results are 
listed below: 

Here the calculated mass is in units of MeV and KgNg means the K and are both 
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TABLE I: Quark shell model calculations 







pure KN 


KN + K*N 


KN + KsNs 


KN + K*N + KsNs 


Sqi 


Parity= — 


2282 


2157 


1943 


1766 


Poi 


Parity= + 


2357.1 


2356.3 


2357.0 


2336.8 



in color octets but coupled to an overall color singlet. The S'-wave state has strong chan- 
nel mixing: the amplitudes of KN, K*N, KsNg, K^Ng are -0.54, -0.29, -0.54, -0.29 
respectively. On the other hand the channel mixing is weak for the P-wave state: The am- 
plitudes oi KN, K*N\, K*N\, K^N^, K^Ng^, K^Ng^ are 0.968, -0.1, ~0, -0.07, 0.224, 
~0 respectively. Here the | and | are the channel spin. The S-wave state definitely has a 
lower mass than that of the P-wave; the channel coupling plays a vital role in reducing the 
calculated S'-wave 0^ mass. In obtaining these results we use the model parameters which 
fit the nucleon mass with five quark components mixed to configuration. Neglecting the 
five quark component, this model will give M{N) = 1.2 GeV, M{K) = 0.8 GeV. Therefore 
it is possible to reduce the penta-quark mass further by taking into account hepta-quark 
component and other channels coupling. We like to point out that our calculation is in fact 
a quark shell model calculation. The large scale shell model calculation of nuclear structure 
shows that it is possible to obtain approximate correct energy and wave function of a nuclear 
state if the Hilbert space is large enough. This should be true for quark shell model as well. 

Lattice QCD and non-perturbative QCD both show that confinement might be due to 
gluon fiux tube (or gluon string) formation in a quark system. The ground state energy of 
the gluon field in a qq meson and q^ baryon can be approximated by a potential j^l 

A 

Vqq = ~ \~ '^qq'^ + Cgg 



— —A^g ^ — — - + Cr^gLmin + C'sg; (2) 

i<j Fi - 

i-'min ^ ^ Li. 

i 

Li is the distance between the quark i and the Y-shaped gluon junction, is the position 
of quark i. The first term in Eq.(0) is the color Coulomb interaction and the second term is 
similar to a linear confinement potential. 

Most of constituent quark models use a quadratic or linear potential to model the quark 
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confinement, 

rij = ri-rj, n = l,2. (3) 

Here A" (a = 1 ■ ■ ■ 8) is the color SU(3) group generator. For a single hadron, qq mesons 
or q^ baryons, such a modelling can be achieved by adjusting the strength constant a of 
the confinement potential. The color factor Aj ■ Xj gives rise to a strength ratio 1/2 for 
baryon and meson which is almost the ratio for the minimum length of the flux tube to the 
circumference of the triangle formed by three valence quarks of a baryon. 

How to extend the confinement potential to multiquark systems is an open question. 



There is a lattice QCD calculation of the penta-quark potential recently j41[. The ground 
state energy of the gluon field in a penta-quark with color structure qq(S)s(3)qq(3) can be 
expressed as 

Vbq = ~ry2 ^ + O'bqLmin + Cso- (4) 

Here qq{3) means a color anti-triplet qq pair. L^m is the minimum length of the color flux 
to connect the five quarks. 

From general SU(3) color group considerations, there might be other color structures for 
a penta-quark: g^(l)gs(l);g^(8)gs(8); gg(3)gg(3)s(3), etc. The first one is a color singlet 
meson-baryon channel; the second is the hidden color meson-baryon channel; the third is 
the color structure used in the Jaffe-Wilczek model. One guesses the energy of these color 
configurations can be expressed in a similar manner as those given in Eq.(2) and (4). Can 
the two body confinement interaction Eq.(3) describe the confinement interaction properly 
for these color structures? It is a question needed to be studied further j^]- A penta-quark 
state should be a mixing of these color structures. Our first model calculation mentioned 
above shows that channel coupling reduces the calculated ground state penta-quark mass. 
It should be true in general that these different color channel mixing will reduce the ground 
state energy. 

To do a model calculation for a multi-quark system with the above multi-body interaction 
and multi-channel coupling is numerically quite involved. We have developed a model, 
called the quark delocalization, color screening model (QDCSM): 1. We re-parameterize the 
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confinement potential Eq.(3) to take into account the effect of cliannels coupling in multi- 
quark systems induced by various color structures, which are not possible for a qq meson and 
q^ baryon; 2. To take into account the orbital excitations but keep the numerical calculation 
simple we still use quark cluster bases but a delocalized quark orbital wave function is used 
to allow the multi-quark system to choose its own favorable configuration, i.e., to allow the 
multi-quark system to vary from the asymptotic hadronic cluster state to a genuine multi- 



quark state and all intermediate configurations |43[. This model explains the existing BB 
interaction data (bound state deuteron and A^A^, A^A, NT, scattering) well with all model 
parameters fixed by hadron spectroscopy except for only one additional parameter, the color 
screening constant /i which is determined by the deuteron properties. More important, it 
is the unique model, so far, which explains a long-standing fact: The nuclear force and the 
molecular force are similar except for the obvious difference of length and ene rgy s cales; the 



nucleus is approximately an A nucleon system rather than a 3A quark system 
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Fig. 2. The effective potentials for 1=0 S-wave. 

The penta-quark mass has been calculated with this model (QDCSM) in a single color 
singlet KN channel approximation. As explained before, the effect of the coupling of other 
color structures and orbital excitations is assumed to have been included in the modelling of 
QDCSM. In the 1 = 1 S-wave KN channel, a pure repulsive effective interaction is obtained 
as shown in Fig.l (corresponding to curve with /i = 1). This is consistent with the KN 
scattering data and helps to rule out the 1 = 1 possibility for the . For comparison 
the naive quark model result is also shown in Fig.l (corresponding to curve with /i = 0), 
which shows a stronger repulsion. In the / = S-wave KN channel, an effective attraction is 
obtained and shown in Fig.2 (/i = 1 curve). This is inconsistent wit 
phase shifts jsJll but might be consistent with Gibbs new results 
MeV is obtained from the minimum of Fig.2 /i = 1 curve. Part of the overestimate of the 
mass is due to the overestimate of K mass, which is 650 MeV in this approach. This can be 
eliminated as follows: One first gets an effective interaction potential from Fig.2 /i = 1 curve 
by subtracting its asymptotic value, then add a zero point oscillating energy r,! (fJ-KN is 
the reduced mass of K and A^, Rq = 0.6/m is the minimum point of the effective potential) 
and the rest mass of and K. In this way one obtains M{<d~^) = 1615MeV . It is still 
about 75 MeV higher than the observed value 1540 MeV. More precise dynamical calculation 
might reduce the mass further. Fig.2 fi = curve shows the naive quark model result 
for comparison. It is almost a pure repulsive interaction and will not accommodate a G"*" 
resonance. The K^n effective potential is shown in Fig.3, which is a very weak repulsive 



1 the VPI KN scattering 
3. A e+ mass of 1706 



interaction resulted from a cancellation of the 7 = 1,0 channel ones. This result shows it is 
hard to get reliable K^n scattering amplitude from K^d scattering data because it is a small 
component in comparison to the big K^p amplitude. Gibbs analysis shows the additional 



complications [l^. To get reliable 1 = scattering phase shifts from K^n is even harder 
because one has to get two big ones (corresponding to / = 1,0 separately), which have 
opposite sign and so cancel each other, from a small one corresponding to (/ = 1) + (/ = 0). 
For the P-wave channels, we only obtain spin averaged effective KN interactions because 
the spin-orbit coupling has not been included yet. In the / = channel, there is a strong 
attraction (shown in Fig.4), as wanted in other quark models with correlations. However 
in our model the P-wave attraction is not strong enough to overcome the kinetic energy 
increase to reduce it to be a ground state. This is consistent with the lattice and QCD sum 
rule results 



36l l38| . In the 7 = 1 channel, only a very weak attraction is obtained. This 



rules out the 1 = 1 0+ again. 
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Fig.4. The effective potentials for 1=0 P-wave. 




Fig.5. The Jaffe-Wilczek configuration. 
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Fig. 6. The energies for Jaffe-Wilczek configuration {/i— 1). 

In a third model, the Jaffe-Wilczek configuration {ud}{ud}s is used but the four non- 
strange quarks are totally antisymmetrized. The space part is fixed to be an equilateral 
triangle with the two diquarks sitting at the bottom corners and the s at top (see Fig. 5). 
The height h and the length L of the bottom side of the triangle are taken as variational 
parameters in addition to the quark delocalization. A three body variational calculation 
with the QDCSM has been done. The minimum of this variational calculation is 1621 MeV 
corresponding to a triangle oi h — 0.6/m, L — 0.6 fm (see Fig. 6). The 0+ mass is simi- 
lar to our second model one. Fig. 7 shows the result obtained with the naive quark model 
Hamiltonian, the minimum is 1799 MeV corresponding to a vanishing triangle. 
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Fig. 7. The energies for Jaffe-Wilczek configuration (/i = 0). 

Based on tliese tliree model results the IJ^ = are both possible to be the ground 
state of the O^. A multi channel coupling quark model calculation with various nontrivial 
color structures are needed to check if the observed 9+ mass can be obtained in quark model 
approach and which parity is the ground state. 

Suppose the is finally verified to be a 1540 MeV narrow width (~ 1 MeV) IJ^ = 0^^ 
state. Then an interesting scenario similar to that of nuclear structure at the 1940's turned 
to the 1950's will recur. The low lying, even parity rotational excitation of nuclei is hard 
to explain by the naive Mayer- Jenson nuclear shell model; Bohr and Mottelson had to 
introduce the rotational excitation of a deformed liquid drop model. Later, nucleon Cooper 
pairs were introduced because of the strong short range pairing correlation. In 1970's-1980's, 
an S-D Cooper pair interaction boson model was developed and the collective rotation was 
re-derived from this model which is based on Mayer- Jenson's nuclear shell model but with 
nucleon pair correlation. In the description of the pentaquark, one has introduced the chiral 
soliton rotational excitation, quark color Cooper pairs and much more. The historical lessons 
of nuclear structure study might be a good pharos to light the way for the study of hadron 
structure. 

Summary: Multi-quark states have been studied for about 30 years. The G"*", if further 
confirmed, will be the first example. Once the multi-quark "Pandora's box" is opened, the 
other multi-quark states: tetra-quark, hexa-quark (or dibaryon), etc., can no longer be kept 
inside. One expects they will be discovered sooner or later and there are claims that some 



tetra-quark states have been observed 



A new landscape of hadron physics 



will appear and it will not only show new forms of hadronic matter but will also exhibit new 
features of low energy QCD. 

Nonperturbative and lattice QCD have revealed the color fiux tube (or string) structure 
of the qq, q^ and even q^q states. The multi-quark system will have more color structures. 
How do these color structures interplay within a multi-quark state? Nuclear structure seems 
to be understood in terms of colorless nucleons within a nucleus. Multi-quark states might 
be not so. We emphasized that the effect of non-trivial color structures of multi-quark 
system should be studied. The low mass and narrow width of the G"*" might be related to 
such new structures instead of to residual interactions. 
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